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Abstract 

This article deals with the problem of optimal allocation of capital to corporate 
bonds in fixed income portfolios when there is the possibility of correlated defaults. 
Under fairly general assumptions for the distribution of the total net assets of a set 
of firms we show that retaining the first few moments of the portfolio default loss 
distribution gives an extremely good approximation to the full solution of the asset 
allocation problem. We provide detailed results on the convergence of the moment 
expansion. We also provide explicit results for the inverse problem, i.e. for a given 
allocation to the set of risky bonds, what is the average risk premium required to make 
the portfolio optimal. Numerous numerical illustrations exhibit the results for realistic 
portfolios and utility functions. 



1 Introduction 

Investors routinely look to the corporate bond market in particular, and spread markets 
in general, to enhance the performance of their portfolios. However, for every source of 
excess return over the risk-free rate, there is a source of excess risk. When sources of 
risk are correlated, the allocation decision to the risky sectors, as well as allocation to 
particular securities in that sector, can be substantially different from the uncorrelated 
case. Since the joint probability distribution of returns of a set of defaultable bonds varies 
with the joint probabilities of default, recovery fraction for each bond, and the number of 
defaultable bonds, a direct approach to the allocation problem that incorporates all these 
factors completely can only be attempted numerically within the context of a default model. 
This approach can have the short-coming of hiding the intuition behind the asset allocation 
process in practice, which leans very heavily on the quantification of the first few moments, 
such as the mean, variance and skewness. In this paper, we will take the practical approach 
of characterizing the portfolio default loss distribution in terms of its moment expansion. 
Focusing on allocation to corporate bonds (the analysis in this paper can be generalized 
to any risky sector which has securities with discrete payoffs), we will answer the following 
questions: 

• In the presence of correlated defaults, how well does retaining the first few moments 
of the portfolio default loss distribution do, for the portfolio allocation problem, as 
compared to the more intensive full numerical solution? 

• For different choices of correlations, probabilities of defaults, and number of bonds in 
the portfolio, how does the optimal allocation to the risky bonds vary? 

• For a target allocation to the risky bond sector, what is the required excess return 
over the actuarially fair value that is a reasonable compensation for the risk when 
default correlations are present? 

In this paper we consider portfolios consisting of risk free assets at a return r and 
corporate bonds with promised return Cj for firm i, and study how correlations between 
defaults affect the optimal allocation to corporate assets in the portfolio. Let hi denote a 
random variable that takes the value 1 if company i defaults in the time horizon At and 
zero otherwise. Pi = E[hi] is the probability of firm i defaulting and dy = Corr[hi,hj] is 
the default correlation between firms i and j. We will quantify the impact of risk associated 
with losses from corporate defaults at some single period time horizon At [] on portfolio 
allocation. We assume that if company i defaults in the time period At a fraction (1 — Ri) 
of the portfolio's original investment, in bonds of company i is lost0. To compensate the 
investor for default risk, the corporate bond's promised return q is greater than that of 
the risk-free assets in the portfolio. The excess return, or spread of investment in bonds 
of firm i can be decomposed into two parts. The first part of the spread arises as the 

1 In this paper, our numerical analysis will assume At to be one year so that default probabilities quoted 
are annualized ones. 

2 Ri is called the recovery fraction. More explicitly, if the original investment in company i is 1 dollar 
then after the time period At the value of the investment is taken to be Ci + Ri dollars if company i defaults 
and 1 + Ci dollars if it survives. Here a is the promised corporate return. 
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actuarially fair value of assuming the risk of default. The second part, which we call fa, is 
the excess risk premium that compensates the corporate bond investor above and beyond 
the probabilistically fair value. In practice, fa can arise due to a number of features not 
directly related to defaults. For example, traders will partition fa into two pieces, one for 
liquidity and one for event risk, fa = + ej. Given that most low grade bonds may have a 
full percentage point arising simply from liquidity premia, fa can fluctuate to large values 
in periods of credit stress. Liquidity ^ is systematic and one would expect it to be roughly 
equal for a similar class of bonds. Event risk contributes to fa due to the firm's specific 
vulnerability to factors that affect it (e.g. negative press). In periods of stress Pi, li and all 
increase simulataneously, and the recovery rate expectations R{ fall, leading to a spike in the 
overall spread. Since these variables can be highly volatile, the reason behind the portfolio 
approach to managing credit is to minimize the impact of non-systematic event risk in the 
portfolio.^] The risk premium itself is not very stable over time. Empirical research shows 
that the excess risk premium might vary from tens of basis points to hundreds of basis 
points. For instance, in the BB asset class, if we assume a recovery rate of 50% and default 
probability of 2%, the actuarially fair value of the spread is 100 basis points. However, it is 
not uncommon to find actual spreads of the BB class to be 300 bp over treasuries [Altman 
(1989)]. The excess 200 bp of risk premium can be decomposed in any combination of 
liquidity premium and event risk premium, and is best left to the judgment of the market. 
When the liquidity premium component is small compared to event risk premium, we would 
expect that portfolio diversification and the methods of this paper are extremely valuable. 

One other factor needs to be kept in mind when comparing historical spreads to current 
levels. In the late eighties and early nineties, the spread was routinely quoted in terms of 
a treasury benchmark curve. However, the market itself has developed to a point where 
spreads are quoted over both the libor swap rate and the treasury rate, and the swap rate 
has gradually substituted the treasury rate as the risk-free benchmark curve. This has two 
impacts. Firstly, since the swap spread (swap rate minus treasury rate for a given maturity) 
in US is positive (and of the order of 50 bp as of this writing), the excess spread needs to 
be computed as a difference to the swap yield curve. Secondly, the swap-spread itself has 
been very volatile during the last few years, which leads to an added source of non default 
related risk in the spread of corporate bonds when computed against the treasury curve. 
Thus, the 200 bp of residual spread is effectively 150 bp over the swap rate in the BB 
example, of which, for lack of better knowledge, equal amounts may be assumed to arise 
from liquidity and event risk premium over the long term. The allocation decision to risky 
bonds strongly depends on the level of risk-aversion in the investor's utility function, and 
the required spread for a given allocation will go up nonlinearly as risk aversion increases. 

The value of the optimal fraction of the portfolio in corporate bonds, here called a op t, 
cannot be determined without knowing the excess risk premium part of the promised corpo- 
rate returns. The portfolio average of the excess risk premium, fa provides the incentive for 
a risk averse investor to choose corporate securities over risk free assets and so the optimal 
fraction of corporates in a portfolio a op t is a function of \x and vanishes at [i = (since there 
is no incentive for a risk averse investor choose the corporate assets when [x = 0). Inverting 
this functional relationship the portfolio average excess risk premium is a function of the 

3 The authors would like to thank David Hinman of PIMCO for enlightening discussions on this topic. 
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optimal fraction of corporate assets, \x = fi(a op t). 

We explore, using utility functions with constant relative risk aversion, the convergence 
of the moment expansion for a opt (n) and /x(a opf ). Our work indicates that, (for a opt less 
than unity) a op t(fJ>) and n(a op t) are usually determined by the mean, variance and skewness 
of the portfolio default loss probability distribution. The sensitivity to higher moments 
increases as a op t does. Some measures of default risk, for example a VAR analysis^, may be 
more sensitive to the tail of the default loss distribution. We also explore how the optimal 
portfolio allocation scales with the number of firms. If defaults are random the optimal 
fraction of corporates increases without bound as the number of firms in the portfolio goes 
to infinity but when there are default correlations a op t goes to a finite value in that limit. 

Historical evidence suggests that on average default correlations increase with the time 
horizon^]. For example, Lucas (1995) estimates that over one year, two year and five year 
time horizons default correlations between Ba rated firms are 2%, 6% and 15% respectively. 
However, the errors in extracting default correlations from historical data are likely to be 
large since defaults are rare. Also these historical analysis neglect firm specific effects 
that may be very important for portfolios weighted towards a particular economic sector. 
Furthermore, in periods of market stress default probabilities and their correlations increase 
[Das, Freed, Geng and Kapadia (2001)] dramatically. 

There are other sources of event risk associated with corporate securities. For example, 
the markets perception of firm i's probability of default Pi could increase over the time 
horizon At resulting in a reduction in the value of its bonds. For a recent discussion 
on portfolio risk due to downgrade fluctuations see Dynkin, Hyman and Konstantinovsky 
(2002). Here we do not address the issue of risk associated with fluctuations in the credit 
spread but rather focus on the risk associated with losses from actual defaults, which for 
the purposes of this paper, are defined as occurring if the assets of a firm fall below its 
liabilities. 

In the next section a simple model for default is introduced. The model assumes that 
the difference between the stochastic assets and liabilities, dj of firm i is an increasing func- 
tion of a normal random variable^] Zi . Default occurs in the time interval At if fluctuates 
below zero. The joint default probabilities are expressed in terms of the correlations and 
volatilities of the risk variables Zj. In section 3 we set up the portfolio problem. Moments 
of the fractional corporate default loss probability distribution are expressed in terms of de- 
fault probabilities and it is shown how these can be used to determine a op t(n) and n(a op t). 
In section 4 the impact of correlations on the portfolio allocation problem is studied using 
sample portfolios where all the firms have the same probabilities of default and the correla- 
tions between firms are all the same. Concluding remarks are given in section 5. Many of 
the mathematical details are relegated to an appendix. 



4 See for example, Jorion (2001). 

5 Zhou (2001) derives an analytic formula for default correlations in a first passage-time default model 
and finds a similar increase with time horizon. 

6 Note that hi itself is not necessarily normal or lognormal. 
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2 A Model For Default 



Assume the ith firm in the portfolio has stochastic assets Yi and debt F^. In the time interval 
At a firm defaults if its assets fall below the value of its debt. Thus, the probability of firm 
% defaulting, in the time interval At, is the probability that eij = Yi — F% < 0. Suppose 



(l; 



Yi- Fi= gi(zi), 



(1) 



where z% is a normal random variable of zero mean and gi is an increasing function of Zi. 
Zi can be thought of as a "production" factor variable upon which the net assets of the 
firm i depend. We impose no further restrictions on the form of the net assets of the firm. 
For example, for some constants hi and qi if gi(zi) = hi + 2i, then is normal and if 
gi(zi) = hi + qi exp(lj) then en is lognormal. The statement that default occurs if net assets 
fall below zero may be restated in terms of the production factor Zi crossing some threshold 
Ti for firm i, i.e. 

Zi <g { - 1] (0) = -Ti. 

The probability distribution for the normal variables Zi 
variances 



(2) 



and their correlation matrix 



E[ziZi], 
E[ziZj] 



is completely specified by their 

(3) 



E[ziZi\E[zjZj] 



(4) 



Without assuming a particular form for gi we cannot deduce the volatilities <Tj and the 
correlations £y from the fluctuations of the fij's. However, we will see later in this section 
that in this model the default probabilities Pj and their correlations dij do determine Tijoi 
and which therefore have a direct connection to the measures that investors use in 
quantifying security risk. 

This is basically the Merton Model[| [Merton (1974)] applied over the time horizon At. 
In this model the joint probability that n companies, which we choose to label 1, ...n, 
default is given by the following integral of default thresholds and correlations: 



(2vr)2ydeT£ 

where the sum goes over i,j = 
tances to default" are 



dx 



l • • • 



dx r 



exp 



(5) 



, n, the scaled default thresholds, or "equivalent dis- 

(6) 



Ti 

Xi = — , 

(Ti 



and is the inverse of the correlation matrix. For n not too large the integrals in 

equation (||) can be done numerically or, since defaults are rare, analytic results can be 
obtained using the leading terms in an asymptotic expansion of the integral. 

7 No time evolution model for the net assets is specified. Allowing default to occur only at the end of the 
time interval the hi can be thought of as the net assets at that time. It is possible to extend the work of 
this paper to a first passage-time model [Black and Cox (1976), Longstaff and Schwartz (1995), Leland and 
Toft (1996), etc.] where default is associated with the first time that Zi crosses a threshold — Ti. We will 
consider the implications of correlated default using a first passage-time model in a further publication. 
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As we mentioned in the introduction, the random variable hi takes the value 1 if firm 
i defaults in the time horizon At and zero otherwise. The default probabilities in equation 
(|5|) are expectations of products of these random variables, 

Ph...i m = E[h h ■ ■ ■ h im ], when i x ^ i 2 ■ ■ ■ / i m . (7) 

Pi is the probability that firm i defaults in the time period At and Pi lr .. ; i m is the joint 
probability that the m-firms ii,...i m default in the time period At. Since hf = hi it 
follows that the correlation of defaults between two different firms (which we choose to 
label 1 and 2) is, 

dj2 = £[n!n 2 ] - £[ni]£[w 2 ] = P12 - P1P2 _ Pi 2 - P1P2 (g) 

12 ^{E[h\] - Elh^Elhl) - E[n 2 ] 2 ) ^1(1-^1)^2(1-^2) ~ \fP~\P~2 

When the correlation £12 and the default probabilities Pi and P2 are small the default 
correlation d\2 is also small. There is a simple situation where default correlations are 
likely to be large. Suppose company 2 is dependent on company 1 so that if company 1 
defaults we know with certainty that company 2 will default. In that case P\i = Pi and 
(neglecting terms suppressed by powers of the default probabilities) the default correlation 
is, 




Note that company 2 can default for reasons unrelated to the health of company 1 and 
equation @ still holds. Suppose company 1 has default probability P\ = 0.001 and company 
2 has default probability P2 = 0.01. Even though both of these are small equation (|9|) implies 
a large default correlation of 0.33. 

In the next section we consider the problem of portfolio allocation for portfolios con- 
sisting of corporate bonds subject to default risk and risk free assets. The implications of 
default risk are addressed using the model discussed in this section. Other sources of risk, 
for example, systematic risk associated with the liquidity part of the excess risk premium, 
are neglected. 



3 The Portfolio Problem 

The total wealth in a portfolio consisting of risk free assets that return r and corporate 
assets that are subject to default risk, after time At, is 

W = W (V - a)(l + r) + a £ f t (l + a - n,-(l - P;))) , (10) 

where, Wq is the initial wealth, a is the fraction of corporate assets in the portfolio, is 
the return on the i'th corporate asset, Ri is the recovery fraction and /, denotes the initial 
fraction of corporate assets in the portfolio that are in firm i (J2i fi = !)■ I n equation ( p^| ) 
the sum over i goes over all N firms in the portfolio. The fractional corporate default loss 
in the portfolio over the time period At is the random variable, 

i=J2fMi-Ri)- (11) 

i 
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The average fractional default loss is 

E[t\ = Y d hPi{l-Ri 



(12) 



and the fluctuations of the fractional loss I about its average value is, 

Si=l-E[t\. 



(13) 



Since the fractional loss I cannot be negative and cannot be greater than one — E[l] < 51 < 
1 — E[l]. The left hand side of this inequality is typically close to zero and the right hand 
side is typically close to one since the expected fractional loss is usually small. The mean of 
51 is zero and the probability distribution for 51 determines the default risk of the portfolio 
associated with fluctuations of the random variables rtj. The moments of this probability 
distribution are, = E[(5l) m ]. Using equation (^) and the property hf = hi the moments 
v (m) can k e expressed in terms of the joint default probabilities. For the first five moments 
the results are given in the appendix. 

We write the promised corporate return (i.e. the expected value of the corporate return) 
for firm i as, 

(H = r + Pi(l - Ri) + im, (14) 
where fii is the excess risk premium and introduce the portfolio average excess risk premium, 



(15) 



If an investor holding the portfolio is risk averse, when /i = the optimal portfolio has only 
risk free assets i.e., a op t = 0. 

In terms of the quantities we have introduced the total wealth after the time period At 
in equation ([!(]) takes the form, 



W = W (l + r + a/j, - aSl 



(16) 



To find out what value of a is optimal a utility function is introduced which characterizes 
the investor's level of risk aversion. Here we use utility functions of the type U^{W) = W" 1 /j 
which have constant relative risk aversion^, 1 — r y. The optimal fraction of corporates, otopt 
maximizes the expected utility of wealth E[U~*(W)]. Expanding the utility of wealth in a 
power series in 51 and taking the expected value gives 



E[U 7 (W)} = (Wj/ 1 )(l + r + atiY 



1 + 



oo 
m=2 y 



T(m — 7) 
— 7)T(m +1) \ 1 + r + a/i 



v 



(m) 



(17) 



where T is the Euler Gamma function. The approximate optimal value of a obtained from 
truncating the sum in equation ([17|) at the m'th moment is denoted by a m . The focus of 
this paper is on portfolios that are not leveraged and have a op t less than unity. We will 
later see that typically, for such portfolios, the a m converge very quickly to a op t so that for 
practical purposes only a few of the moments need be calculated. 



3 See, for example, Ingersoll (1987). 
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The utility of wealth has explicit dependence on [i so the optimal value of the fraction 
of corporates a op t is a function of it i.e., a op t = a op t(fi). Expanding in a power series about 
/i = 0, 



oo 

cMaO = E -^ n - ( 18 ) 

n=l 



The coefficients Si can be expressed in terms of the moments, v ^ . Explicitly, for the first 
two coefficients^, 

and 

(2 - 7 )(i + r y3) 

32 -~ (l- 7 )2(t,(2))3 • ( 20 ) 

a ^ - = (i- 7 y2) » ( 21 ) 

has the proper behavior as the various parameters it depends on vary. Increasing the 
portfolio average excess risk premium [i increases the fraction of corporates. Decreasing 
7 corresponds to greater risk aversion and gives a lower value of a op t- If the excess risk 
premium [i vanishes then the optimal portfolio contains only risk free assets. 

Inverting equation (|l8|) gives the risk premium as a function of the optimal fraction of 
corporates, fx = n(a op t). Truncating the sum in equation (|lg|) at the second term 



The approximation, 



. »1 - V S 1 + 2a o P tS2 If s 2 2 \ 

V~i{a pt) ^ = = — a opt - ^a opt H I . (22) 



-s 2 si V 2sf 

We have added a subscript 7 to emphasize that the excess risk premium depends on the 
assumed level of risk aversion. Note that equation (^2|) is not a consequence of the moment 
expansion but rather relies on an expansion of the utility of wealth in the excess risk 
premium. If the excess risk premium is large even when the moment expansion for a op t 



works quite well equation (22) may not be useful. 



4 Sample Portfolios 

Here we consider sample portfolios where the correlations and volatilities of the fj and the 
default thresholds are the same for all N firms i.e., £y = £ and \i = X- Then all the 
probabilities of default are the same Pi = P and the joint default probabilities are also 
independent of which firms are being considered, P^...^ = P\2...m- The corporate returns 
are also assumed to be independent of firm i.e., Cj = c and m = fi. We also take the 
portfolios to contain equal assets in the firms so that, /, = 1/N, for all i, and assume all 
the recovery fractions are zero. 

9 A derivation of these results for si and S2 is given in the appendix. 
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There are two cases where results are easily available without using the moment expan- 
sion of the utility function. The first of these is random defaults where the correlations van- 
ish. In this case the probability of a fraction loss of I = n/N is (1 — P) N ~ n P n N\ / (N — n)\n\ 
and so the m'th moment of the loss distribution is 

^ Q \N J (N-n)\n\ 
The expected utility of wealth is 

E[U,{W)\ = «/ 7 ) - P) N ~ n P n {N m n)]nl (l + r + an - a (± - P )) 7 . (24) 

Having the explicit expression for the utility of wealth in equation ( p4|) lets us compare 
results of the moment expansion for the optimal fraction of corporates a m with the all 
orders result, a opt - These are presented in Table I in the case where the probability of 
default is P = 2%, the risk free return is 3.5% . According to equation ( |l4[ ) this implies a 
total promised corporate return of c = 5.6%. Results for different values of the number of 
firms N and the level of risk aversion 1 — 7 are shown in Table I. We also give in columns 
three and four of Table I the volatility, vol = V v^ 2 \ and the skewness, skew = /Vv^ , 
of the portfolio default loss distribution. 

Increasing N gives a larger value of the optimal fraction of corporates because diversi- 
fication reduces risk. This occurs very rapidly with N. For N = 1, fj, = lOObp and 7 = —4 
the optimal fraction of corporates is only 8.3%. By N = 10 increasing the number of firms 
has reduced the portfolio default risk so much that the lOObp excess risk premium causes a 
portfolio that is 83% corporates to be preferred. 

For all the entries in Table I the moment expansion converges very rapidly, although for 
low iV it is the small value of a op t that is driving the convergence. Since in equation (|17|) the 
term proportional to t/ m * > has a factor of a m accompanying it we expect good convergence 
of the moment expansion at small a. 

Table I: Optimal Fraction of Corporates for Sample Portfolios with Random (£ = 0) 
Defaults and Probability of Default P = 2% 



7 


N 


vol x 10 3 


skew 


/' 


a 2 


«3 


«4 




a opt 


/0.01 


1 


140 


6.9 


lObp 


0.053 


0.051 


0.051 


0.051 


0.051 \ 


0.01 


5 


63 


3.1 


lObp 


0.27 


0.25 


0.25 


0.25 


0.25 


0.01 


10 


44 


2.2 


lObp 


0.53 


0.51 


0.51 


0.51 


0.51 


0.01 


15 


36 


1.8 


lObp 


0.80 


0.77 


0.76 


0.76 


0.76 


-4 


1 


140 


6.9 


lOObp 


0.11 


0.086 


0.083 


0.083 


0.083 


-4 


5 


63 


3.1 


lOObp 


0.54 


0.43 


0.42 


0.41 


0.41 


I" 4 


10 


44 


2.2 


lOObp 


1.1 


0.88 


0.84 


0.83 


0.83 j 



The focus of this paper is on portfolios that are not leveraged and have a op t < 1 . But a 
value a op t > 1 is not forbidden when finding the maximum of the expected utility of wealth. 
This occurs at lowest order in the moment expansion in the last row of Table I. 
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The other situation where the optimal portfolio can be found without using the moment 
expansion is when the correlations are maximal, i.e. £ = 1. Then all the companies must 
default together and so for any N the utility of wealth is given by equation (^4|) with N = 1. 
Explicitly, for maximal correlations, 

E[U 7 (W)] = (W 7 / 7 ) [(1 - P) (1 + r + an + aPf 

+ P {I + r + an - a{l - P)) 1 } . (25) 

Correlations make diversification less effective and in the extreme case of maximal correla- 
tions the optimal fraction of corporates is independent of the number of firms. 

Including correlations makes the moment expansion converge slower. The situation 
will never be worse than the maximal correlation case and so it is worth examining its 
convergence. This is done in Table II. Here we take, £ = 1, r = 0.035 and P = 2.0%. 
Results are presented for 7 = 0.01 and 7 = —4 and excess risk premiums between lObp and 
400bp. 

With 7 = 0.01 the excess risk premium ju(54%) is 200bp and for that case 05 is within 
4% of the correct value. Just including the variance v ^ gives you a result that is off by 
almost a factor of two but 03 is within 30% of the correct value. In realistic cases where 
correlations are significantly less than unity the convergence of the moment expansion should 
be substantially better. 

Note that for 7 = 0.01 and excess risk premium of 400bp, gives an optimal portfolio 
with almost 75% corporates. This excess risk premium is too small for most investors to 
risk, with a 2% probability, losing 75% of their wealth. With 7 = —4 the investor requires 
an excess risk premium of 46.6% before the optimal fraction of corporates reaches 75%. 

Table II: Worst Case Convergence of Moment Expansion, £ = 1, for Probability of Default 

P = 2%. 



7 


vol x 10 3 


skew 


A* 




«3 


OCA 


«5 


C^opt 


/0.01 


140 


6.9 


lObp 


0.053 


0.051 


0.051 


0.051 


0.051 \ 


0.01 


140 


6.9 


lOObp 


0.54 


0.39 


0.37 


0.36 


0.36 


0.01 


140 


6.9 


200bp 


1.1 


0.68 


0.59 


0.56 


0.54 


0.01 


140 


6.9 


300bp 


1.8 


0.92 


0.77 


0.70 


0.66 


0.01 


140 


6.9 


400bp 


2.6 


1.2 


0.92 


0.83 


0.74 


-4 


140 


6.9 


lOObp 


0.11 


0.086 


0.083 


0.083 


0.083 


-4 


140 


6.9 


200bp 


0.21 


0.15 


0.14 


0.14 


0.14 


-4 


140 


6.9 


300bp 


0.33 


0.21 


0.19 


0.18 


0.18 


V "4 


140 


6.9 


400bp 


0.45 


0.26 


0.23 


0.22 


0.21 / 



Finally we consider the more typical case £ ^ 0, 1. The formulae used to determine the 
joint default probabilities and the moments — are given in the appendix and using 
them we calculate 02 — 05. 
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Table Ilia: Moment Expansion for Optimal Portfolio with Correlated Defaults. The 

Probability of Default is P = 2%. 



7 


AT 


K 




VOl X 1U 


skew 




Ol2 


«3 


04 


05 




/0.01 


10 


0.50 


0.152 


68 


5.1 


lObp 


0.22 


0.21 


0.21 


0.21 




0.01 


10 


0.45 


0.126 


64 


4.8 


lObp 


0.25 


0.23 


0.23 


0.23 




0.01 


10 


0.40 


0.103 


61 


4.5 


lObp 


0.28 


0.26 


0.26 


0.26 




0.01 


10 


0.35 


0.0823 


58 


4.2 


lObp 


0.31 


0.29 


0.29 


0.29 




0.01 


10 


0.30 


0.0645 


56 


3.8 


lObp 


0.34 


0.32 


0.32 


0.32 




0.01 


10 


0.25 


0.0491 


53 


3.5 


lObp 


0.37 


0.35 


0.35 


0.35 




0.01 


50 


0.50 


0.152 


58 


5.5 


lObp 


0.32 


0.29 


0.29 


0.29 




0.01 


50 


0.45 


0.126 


53 


5.2 


lObp 


0.37 


0.34 


0.34 


0.34 




0.01 


50 


0.40 


0.103 


49 


4.8 


lObp 


0.44 


0.41 


0.40 


0.40 




0.01 


50 


0.35 


0.0823 


44 


4.4 


lObp 


0.53 


0.49 


0.48 


0.48 




0.01 


50 


0.30 


0.0645 


40 


4.0 


lObp 


0.64 


0.59 


0.58 


0.58 




0.01 


50 


0.25 


0.0491 


37 


3.5 


lObp 


0.37 


0.78 


0.72 


0.71 




0.01 


100 


0.50 


0.152 


56 


5.6 


lObp 


0.33 


0.30 


0.30 


0.30 




0.01 


100 


0.45 


0.126 


51 


5.3 


lObp 


0.40 


0.36 


0.36 


0.36 




0.01 


100 


0.40 


0.103 


47 


4.9 


lObp 


0.48 


0.44 


0.43 


0.43 




0.01 


100 


0.35 


0.0823 


42 


4.6 


lObp 


0.58 


0.53 


0.52 


0.52 




0.01 


100 


0.30 


0.0645 


38 


4.1 


lObp 


0.72 


0.66 


0.65 


0.65 




Vo.oi 


100 


0.25 


0.0491 


34 


3.6 


lObp 


0.91 


0.83 


0.82 


0.82 


/ 



Table Illb: Moment Expansion for Optimal Portfolio with Correlated Defaults. The 

Probability of Default is P = 2%. 





7 


N 






vol x 10 3 


skew 




«2 


03 


04 


«5 




/ 


-4 


10 


0.50 


0.152 


68 


5.1 


lOObp 


0.45 


0.34 


0.31 


0.31 






-4 


10 


0.45 


0.126 


64 


4.8 


lOObp 


0.50 


0.38 


0.35 


0.34 






-4 


10 


0.40 


0.103 


61 


4.5 


lOObp 


0.56 


0.42 


0.39 


0.38 






-4 


10 


0.35 


0.0823 


58 


4.2 


lOObp 


0.62 


0.47 


0.44 


0.43 






-4 


10 


0.30 


0.0645 


56 


3.8 


lOObp 


0.69 


0.52 


0.48 


0.47 






-4 


10 


0.25 


0.0491 


53 


3.5 


lOObp 


0.75 


0.57 


0.54 


0.53 






-4 


50 


0.50 


0.152 


58 


5.5 


lOObp 


0.64 


0.45 


0.41 


0.40 






-4 


50 


0.45 


0.126 


53 


5.2 


lOObp 


0.76 


0.53 


0.48 


0.47 






-4 


50 


0.40 


0.103 


49 


4.8 


lOObp 


0.91 


0.63 


0.57 


0.55 






-4 


50 


0.35 


0.0823 


44 


4.4 


lOObp 


1.1 


0.76 


0.68 


0.66 






-4 


50 


0.30 


0.0645 


40 


4.0 


lOObp 


1.3 


0.92 


0.83 


0.79 






-4 


50 


0.25 


0.0491 


37 


3.5 


lOObp 


1.7 


1.1 


1.0 


0.98 






-4 


100 


0.50 


0.152 


56 


5.6 


lOObp 


0.67 


0.47 


0.43 


0.41 






-4 


100 


0.45 


0.126 


51 


5.3 


lOObp 


0.81 


0.56 


0.51 


0.49 






-4 


100 


0.40 


0.103 


47 


4.9 


lOObp 


0.98 


0.67 


0.60 


0.68 






-4 


100 


0.35 


0.0823 


42 


4.6 


lOObp 


1.2 


0.99 


0.95 


0.94 






-4 


100 


0.30 


0.0645 


38 


4.1 


lOObp 


1.5 


1.0 


0.91 


0.86 




V 


-4 


100 


0.25 


0.0491 


34 


3.6 


lOObp 


1.9 


1.3 


1.1 


1.1 


/ 
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Tables III give the results of this calculation for portfolios with r = 3.5% and P = 2%. 
The levels of risk aversion used are 7 = 0.01 and 7 = —4 and respective excess risk premiums 
are taken to be \i = lObp and fj, = lOObp. The number of firms in the portfolios are taken to 
be N = 10, 50 and 100. Unlike Table I the defaults are not random and values of «2 — «5 axe 
given for correlations £ ranging from 0.50 to 0.25. The default correlations d y - corresponding 
to these choices of £ are listed in the fourth column of Tables III. The volatility and skewness 
of the portfolio default loss distribution are given in columns five and six. The joint default 
probabilities used to construct Tables III are given in the appendix. 

For all the entries in Tables III the convergence of the moment expansion is good. In 
most cases the second moment gives a reasonable approximation and by the third moment 
it is usually quite accurate. The convergence is better at lower risk aversion and lower 
values of a op t- 

Correlations dramatically effect the dependence of the optimal portfolio allocation on the 
total number of firms. For the £ = 0.5 entries in Table Ilia the optimal fraction of corporates 
is, 0.21, 0.29, and 0.30 for N = 10, 50 and 100 respectively. Increasing the number of firms 
beyond 50 only results in a small increase in the optimal fraction of corporates. For random 
defaults (i.e., £ = 0) the optimal fraction of corporates goes to infinity as the number of 
firms in the portfolio goes to infinity. That is because the moments of the portfolio loss 
distribution go to zero as N — > 00. For example, with £ = the variance of the default loss 
distribution is, 

J* = (26) 
and the skewness of the default loss distribution is 

v ® 1 

[1-2P). (27) 



(„(2))3/2 ^NP(1 - P) 



For random defaults as N — > 00 the distribution for 81 approaches the trivial one where 
81 = occurs with unit probability. However for £ > the moments v^ m ' go to non- 
zero values in the limit N — > 00 and the default loss distribution remains non-trivial and 
non-normal. 

The examples in this section have all used a default probability of 2%, however, the 
convergence of the moment expansion is similar for significantly larger default probabilities. 
Using r = 3.5%, P = 10%, 7 = -4, N = 100 \i = 300bp and £ = 0.25 we find that: 
a<i = 0.66, 03 = 0.51, 04 = 0.48 and a*, = 0.47 Again the convergence of the moment 
expansion is quite good. 

The impact of changing the probability of default on the portfolio average excess risk 
premium needed to have an optimal portfolio with 25% corporates is studied in Table IV. 
It contains approximate values of fi(25%) for default probabilities P = 2%, 4%, 6%, 8% 
and 10%. Other parameters are: r = 3.5%, £ = 0.25 and 7 = -4. In Table IV fi^(25%) 
is the excess risk premium that follows from just keeping the variance (i.e. only s\ in 
equation (p2|)) and fi^(25%) is the value that follows from keeping both the variance and 



skewness (i.e. both s\ and S2 in equation (22)). Note that even though we are holding the 



'production" factor variable correlations £ fixed in Table IV the default correlations change 
because we are changing the probabilities of default. 
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Correlations change the behavior of n(a op t) as the number of firms in the portfolio gets 
large. For uncorrelated defaults the portfolio average excess risk premium goes to zero as 
N — > oo but with correlations it goes to a nonzero value in this limit. For the P = 10% 
entries in Table IV the required excess risk premium to achieve an optimal portfolio with 
25% corporates drops from 264bp to 157bp to 145bp in going from 10 to 50 to 100 firms. For 
N = 10, 000 it is 133bp indicating that going beyond 100 firms can further reduce /i(25%) 
by only about 10%. 



Table IV: Approximate Values for //(25%) for Different Default Probabilities. 



7 


N 


p 


£ 


dij 


vol x 10 3 


skew 


/i«(25%) 


/x( 2 )(25%) 


-4 


10 


2% 


0.25 


0.0491 


53 


3.5 


34bp 


41bp 


-4 


10 


4% 


0.25 


0.0691 


79 


2.6 


75bp 


92bp 


-4 


10 


6% 


0.25 


0.0833 


99 


2.2 


120bp 


149bp 


-4 


10 


8% 


0.25 


0.0945 


117 


1.9 


164bp 


207bp 


-4 


10 


10% 


0.25 


0.104 


132 


1.7 


210bp 


264bp 


-4 


50 


2% 


0.25 


0.0491 


37 


3.5 


16bp 


18bp 


-4 


50 


4% 


0.25 


0.0691 


58 


2.7 


41bp 


47bp 


-4 


50 


6% 


0.25 


0.0833 


76 


2.3 


69bp 


81bp 


-4 


50 


8% 


0.25 


0.945 


91 


2.0 


lOObp 


118bp 


-4 


50 


10% 


0.25 


0.104 


105 


1.8 


132bp 


157bp 


-4 


100 


2% 


0.25 


0.0491 


34 


3.6 


15bp 


5.9bp 


-4 


100 


4% 


0.25 


0.0691 


55 


2.8 


36bp 


42bp 


-4 


100 


6% 


0.25 


0.0833 


72 


2.3 


63bp 


73bp 


-4 


100 


8% 


0.25 


0.0945 


87 


2.0 


92bp 


108bp 


-4 


100 


10% 


0.25 


0.104 


101 


1.8 


122bp 


145bp 



The approximate values of n(a op t) in Table IV are close to the correct values. For example, 
taking P = 10% and = 100 we find that an excess risk premium [i = 145bp gives: 
a 2 = 0.30, a 3 = 0.26, a 4 = 0.26 and a 5 = 0.26. 

For 7 = —4, a op t > 25% and £ = 0.25, an expansion in the risk premium does not 
converge fast enough for equation (^) to be useful for /J,(a op t). However, one can still find 
fJ>(otopt) using the moment expansion for a op t- For example with 7 = —4, A?" = 100, £ = 0.25, 
r = 3.5% and P = 10% we find that /x(50%) = 330bp. 

5 Concluding Remarks 

Default correlations have an important impact on portfolio default risk. Given a default 
model the tail of the default loss probability distribution is difficult to compute, often in- 
volving numerical simulation) 10 ! of rare events. The first few moments of the default loss 
distribution are easier to calculate, and in the default model we used this involved some 
simple numerical integration. More significantly, the first few terms in the moment expan- 
sion have a familiar meaning. Investors are used to working with the classic mean, variance 

10 See, for example, Duffie and Singleton (1999) and Das, Fong and Geng (2001). 
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and skewness measures and have developed intuition for them and confidence in them. In 
this paper we studied the utility of the first few moments of the default loss probability 
distribution for the portfolio allocation problem. 

The default model we use assumes that the assets minus liabilities of firm i, hi, is 
a function of a normal random variable. However, the cij themselves are not necessarily 
normal or lognormal and their probability distribution can have fat tails. When the fij are 
correlated the portfolio default loss distribution does not approach a trivial)^] or normal 
probability distribution as the number of firms in the portfolio goes to infinity. Correlations 
dramatically decrease the effectiveness of increasing the number of firms, N, in the portfolio 
to reduce portfolio default risk. 

The value of the optimal fraction of corporate assets, a op t, cannot be determined without 
knowing the portfolio average of the excess risk premium part of the corporate returns, fi. 
It provides the incentive for a risk averse investor to choose corporate assets over risk free 
assets and so a op t is a function of \x which vanishes at [i = 0. Inverting this functional 
relationship the portfolio average excess risk premium is a function of the optimal fraction 
of corporate assets, (i = fi(a op t). This is the amount of "free money" over the actuarially 
fair value of the investment an investor that desires an optimal portfolio requires before 
taking on the added risk of having a fraction a op t of corporate securities in the portfolio. 
We explored the convergence of the moment expansion for a op t(/j,) and /i(a op t). Our work 
indicates that, for a op t less than unity, the convergence of the moment expansion is quite 
good. The convergence of the moment expansion gets poorer as a op t gets larger and as the 
level of risk aversion gets larger. 

The values of a op t(/i) and n(a opt ) depend on the utility function used. In this paper we 
used utility functions with constant relative risk aversion, 1 — 7. It is possible to make other 
choices and while the quantitative results will be different for most practical purposes we 
expect that the general qualitative results should continue to hold. 
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A Selected Mathematical Details 

The moments of probability distribution for fluctuations in the fractional default loss are 
defined by, 

„M = E[(5l) m }. (28) 
It is convenient to introduce the quantity, 

fi = Ml - Ri). (29) 

Explicit expressions for the first five moments in terms of the default probabilities, Pi lt ...i m 
are, 

v (2) = E fifjPii + E fl p i - E $\ 2 > ( 3 °) 

v {3) = E fifJkPijk + 3 E fifj p ij + E fi p i ~ ^E\l] - E\l]\ (31) 

v (A) = E fifjfkliPijki + e E ~fUJkPijk+^~ffh p i3 

+ 3^fifJP ij + Y:fiPi-^E\l}-6v^E\l} 2 -E\l}\ (32) 

i^j i 
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and 

v {5) _ fifjfkflfmPijklm + 10 ^2 fi fjfkflPijkl + 10 fifjfkPijk 

+ 15 Yl fifjfkPijk + io J2fifJ p v + 5 Y,fifj p ij + Y,fi p i 

i++3+k i+j i^j i 

- 5v^E[l] - 10v<® E[t\ 2 - Wv^E[l] 3 - E[l] 5 . (33) 

Even if the defaults are random the probability distribution for 81 has non zero higher 
moments. As the number of firms in the portfolio increases the importance of this goes 
down in comparison with the effects of correlations. 

The expression for the variance, v^ 2 \ given there can be rewritten as, 



#3 



>'- -EM(^-^) + EM-^ (34) 



Equation (|34j ) can be expressed in terms of the default probabilities Pi and default correla- 
tions dij. Assuming that the Pi are small, 



„(2) ~ 



Y,fifj\l p i p j d ij+Y,~ft p A- ( 35 ) 



The maximum of the utility of wealth occurs at a value of a where its first derivative 
vanishes. Differentiating the utility of wealth with respect to a and equating it to zero 
gives, 

„ _M_ _ (1~7V 2 W _ (2-7)(l-7)tt (3) <f (m 
1 + r (1+r) 2 2(1 + r) 3 "" 1 ' 



Using the expansion of a op t is powers of the excess risk premium given in equation ([18]) this 
becomes, 

(1 - 7)t/ 2 )si / u (1 - 7)v( 2 )s 2 ^ 2 



1 + r (1+r) 2 2(1 +r) 

(2-7)(l -7)f (3) sf/x 2 



2(1 + r) ; 



+ ..., (37) 



where the ellipses denote terms higher orderP|in the excess risk premium [i. To derive 
equations di~9| ) and (^) we note that since fi is arbitrary the coefficients of \i and fj? must 
vanish yielding, 

51 = (1-7)^) ' (38) 

and 

_ (2 - 7 )s 2 vW __ (2-7)(l + r)t/ 3 ) 

52 ?;( 2 )( l + r) (1- 7 )V 2 )) 3 ' 1 ' 

12 In equation ( p6[ ) the ellipses denote terms that give rise to effects higher than order /j 2 in equation (p7[). 
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For sample portfolios discussed in section 4 all the firms have the same probability of 
default and fi = fi = 1/N. In that case, to express the moments v^ m ' in terms of the 
default probabilities we use equations (p0|) to (5^) which yield the following results, 



„(2, = (iV^l) pi2 + ^ p _ p2j {40) 
^ _ W-W-*) ^ + 3 (^i) pi2 + l_ ip _ 3 „ (2)p _ p 3 (41) 



(4) (jy-i)(JV-2)(JV-3) „ (Ar_i )(A r_ 2) (jy^) 

1> U = ^3 P1234 + 6 ^ Pl23 + 7^^Pl 2 

+ -^P-4v^P-6vWp 2 -P 4 , (42) 



and 



(5) (jV-l)(jV-2)(jV-3)(jV-4) n (iV -l)(JV-2)(JV-3) „ 

v = at4 -^12345 + 1U ^ Pl234 

(N-l)(N-2) (N-l)P 12 1 a -\ 

- 10w (3) P 2 - 10?/ 2 )p 3 - P 5 . (43) 

The joint default probabilities are given by the integrals in equation (|5|) which we eval- 
uate by numerical integration. For this we need the inverse and determinant of a m x m 
correlation matrix (m < N) with off diagonal elements £. It has one eigenvalue^ equal to 
1 + (m — 1)£ and the others equal to 1 — £. Consequently its determinant is 

det[^] = (l-O m - 1 [l + (^-lK]. (44) 
Its inverse has diagonal elements, 



,(_!) _ 1 + (m - 2)g 

1 + (m - 2)£ - (m - 1)£ 2 



£7^ = i^A^ tTi (45) 



and off diagonal elements (i 7^ j), 



£ { ~ 1} = (46) 

?4J 1 + (m - 2)f - (m - 1)£ 2 k j 

In Table V the joint default probabilities are given for P = 2%. 

13 For large portfolios this correlation matrix is only consistent for £ non-negative since a correlation 
matrix must have non-negative eigenvalues. It is very important that the correlation matrix £y is positive 
semi-definite. If it has negative eigenvalues the integrals in equations (^) are not well defined. Typically 
a correlation matrix that is forecast using qualitative methods will have some negative eigenvalues. A 
practical method for constructing the consistent correlation matrix that is closest to a forecasted one is 
given in Rebonato and Jackel (2000). 
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Table V: Joint Default Probabilities when P = 2%. 
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K A 1 

0.41 
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z.o4 


1.51 


i no 
1.02 


0.O5 


/I ^ A 

4.04 


i.yo 


l.Ui 


U.O (\) 


U.OU 


3.98 


1.49 


0.746 


0.439 


0.55 


3.39 


1.12 


0.504 


0.273 


0.50 


2.87 


0.820 


0.330 


0.162 


0.45 


2.41 


0.589 


0.208 


9.13 x 10~ 2 


0.40 


2.01 


0.412 


0.124 


4.79 x 10~ 2 


0.35 


1.66 


0.279 


7.01 x 10~ 2 


2.30 x lO" 2 


0.30 


1.36 


0.181 


3.67 x 10^2 


9.85 x 10~ 3 


0.25/ 



Using equations p0| ) to and the results of Table V the moments of the default loss 
probability distribution can be calculated. They are given in Table VI for the case N = 100. 

Table VI: Moments of the Default Loss Probability Distribution when P = 2% and 

N = 100. 



(2) x 10 3 


v (3) x 10 3 


t/ 4 ) x 10 3 


u(5) X 10 3 




11.3 


8.14 


6.66 


5.71 


0.90 \ 


9.55 


6.29 


4.85 


3.97 


0.85 


8.17 


4.91 


3.56 


2.79 


0.80 


7.01 


3.84 


2.62 


1.95 


0.75 


6.01 


2.99 


1.91 


1.35 


0.70 


5.15 


2.31 


1.38 


0.923 


0.65 


4.40 


1.77 


0.975 


0.614 


0.60 


3.74 


1.33 


0.675 


0.397 


0.55 


3.15 


0.988 


0.455 


0.247 


0.50 


2.64 


0.717 


0.297 


0.147 


0.45 


2.19 


0.506 


0.185 


8.27 x 10~ 2 


0.40 


1.79 


0.346 


0.110 


4.34 x 10~ 2 


0.35 


1.45 


0.227 


6.16 x 10~ 2 


2.09 x 10~ 2 


0.30 


1.15 


0.142 


3.18 x 10~ 2 


8.97 x 10~ 3 


0.25 j 



For £ = 0.25 the skewness of the default loss probability distribution is 3.6 indicating that 
for £ > 0.25 this probability distribution is highly nontrivial and is very far from a normal 
distribution. Despite this the moment expansion is & useful tool for determining ot pt. 
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